Monopole and multipole plasmons in a two-dimensional system 
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Using monolayer graphene as a model system for a purely two-dimensional (2D) electron gas, 
we show by energy electron loss spectroscopy, highly resolved both in energy and momentum, that 
there is a significant probability for the excitation of not only one but two dispersing losses. The 
appearance of both losses is independent of the substrate (we tested graphene on the Si face of 6H- 
SiC(OOOl), and on Ir(lll) without and with intercalated Na layer), and the ratio of the slope in the 
dispersion curves varies between 1.4 (SiC) and 2. While the lower dispersion curve can be attributed 
to the excitation of the monopole plasmon, in agreement with theoretical model calculations, the 
upper dispersion branch has not been identified before for plasmonic excitations in a 2D electron 
gas, and we assign it to the excitation of a multipole sheet plasmon. 

PACS numbers: 73.22.Lp 73.22.Pr 79.20.Uv 



I. INTRODUCTION 

The properties of sheet or wire plasmons in low- 
dimensional systema^ differ strongly from those in three 
dimensions. They are characterized by an acoustic-like 
dispersion, i.e., their excitation energy goes to zero in 
the limit of large wavelength^^lSl_ This is due to the lack 
of a restoring force for charge density oscillations in the 
long wavelength limit, contrary to the behavior of surface 
plasmons on the surfaces of bulk metals, which originate 
simply from the truncation of the bulk electronic system. 

These surface plasmons have been identified as self- 
sustained collective excitations at metal surfaced long 
ago, but they gained renewed interest in recent years by 
investigations of the electromagnetic properties of nanos- 
tructured materiald^ and perspectives for many poten- 
tial applications, e.g. in conjunction with engineering 
cavity-free and broad band photon-emitter interactionglS' 
via subwavelength co nfinem ent of optical fields near 
metallic nanostructureJ^^'^. The use of sheet plasmons 
in this context, because of their much flatter dispersion 
compared with "conventional" surface plasmons, would 
allow an ex- treme confinement of light in a broad fre- 
quency range from medium to far infrared. 

In this paper we want to put emphasis on the fact that 
there is more than one dispersing plasmon mode even in 
the plasmons of 2D sheets. This is well known for metal- 
lic surfaces of 3D systems, but such a mode has not been 
identified for low-dimensional systems yet. For metal- 
lic bulk systems, an additional plasmonic surface mode 
has first been found theoretically assuming a smoothly 
decreasing electron density profile through the surface^^, 
which appeared at higher energies than the standard lon- 
gitudinal mode. This is the so-called multipole surface 
plasmon, which was experimentally corroborated only 
much latei^i^lti^, mostly on the surfaces of alkali metals 
and on Al, and by quantitative DFT calculation^i^^Hi^. It 
should be pointed out here that the intensity and even the 
detectability of the multipole plasmon depends explicitly 
on the shape of the electronic density distribution normal 
to the surfac^^ni. 



The adsorption of a 2D conducting sheet, e.g. in form 
of a graphene monolayer on SiC(OOOl) or on Ir(lll), rep- 
resents a localized electron distribution, whose density 
decays exponentially to both sides normal to the surface. 
This decay will be asymmetric because of the screening 
properties of the substrate, described by its dielectric 
function. Therefore, they are also candidates for the ex- 
istence of multipole plasmons, but surprisingly, they have 
not been identified so far. Only asymmetric peak shapes 
have been detected in most cases'^^"^^. With simultane- 
ous high energy and momentum resolution, however, the 
multipole plasmon loss can be identified, as we will show 
below. 

Due to the low conductance of Id and 2d metallic sys- 
tems, one has to cope with comparatively strong damp- 
ing of plasmonic excitations and their sensitivity to de- 
fects on the atomic scale. Therefore, such studies re- 
quire geometrically perfect and at the same time elec- 
tronically flexible materials. Graphene with its two- 
dimensional lattice of sp^-hybridized carbon atoms has 
been shown to be an almost ideal model system for 
studies of many fundamental aspects of 2d electron sys- 
tems, e.g. their electronic correlations, collective phe- 
nomena, many-body interactions, dynamical processes 
and their interrelation^2lH2Zl_ Single layers can be synthe- 
sized by different techniques, e.g. exfoliation of graphite 
or decomposition of hydrocarbons on transition metal 
surfaces^"' SiC substrates by sublimation of Si allow 

detailed studies of the morphology, the interface, the elec- 
tronic structure, and even of transport propertied^^HMI, 
The decay mechanisms of plasmons in graphene layers, 
the consenquences on the form of the dispersion curve, 
and their sensitivity to imperfections has been subject 
of several studies by our groufPlHlZl While these stud- 
ies show that the decay mechanisms have to be explicitly 
taken into account for the understanding of the disper- 
sion of the plasmons in 2D systems like graphene due to 
their short plasmonic life times, wc will concentrate here 
on a detailed line shape analysis of the loss peaks and a 
comparison of the results obtained with single graphene 
sheets on SiC(OOOl) and on Ir(lll) both on the clean 
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substrate and intercalated with a submonolayer of Na. 
This comparison will enable us to see also differences by 
the changing position of the Fermi level in graphene and 
get new insights into the excitation mechanisms of low- 
dimensional plasmons. 



II. EXPERIMENT 

We performed all measurements under UHV conditions 
at a base pressure of 5 x 10~^Pa. As substrates, Si- 
terminated 6H-SiC(0001) samples (n-doped, ~ lO^^cm"'^ 
from SiCrystal AG) were used, which were etched by an 
ex-situ procedure in a furnace of H2 atmosphere in order 
to remove the residual roughness from polishing steps 
[29] . Graphene was grown by sublimation of Si while an- 
nealing the SiC substrate to approximately 1500 K and 
was stopped after generation of the buffer layer and one 
additional graphene layer. Further details can be found 
i] j35 | 36 | These procedures resulted large flat terraces with 
an average terrace width of about 150 nm, as checked 
with high resolution LEEE^^. 

For graphene on Ir(lll) clean and well ordered Ir sur- 
faces were generated in a different vacuum system'^^ by 
repeated cycles of Kr+-sputtering at 1100 K followed by 
subsequent flash annealing to 1500 K. Graphene mono- 
layers were grown on these surfaces by decomposition 
of hydrocarbons (ethylene) at temperatures around 1320 
K, at a pressure of 2 x 10"'' Pa for 20 sec. More de- 
tails are described elsewher^^S, These samples were then 
transferred into the ELS-LEED systems. Mild annealing 
to 500 K for several minutes leads to full desorption of 
water and other contaminants, as obvious from a perfect 
diffraction pattern of a graphene monolayer structure ob- 
tained with SPALEED. The spot profile analysis of the 
(00) beam in LEED revealed the terrace width on this 
surface was approximately 75 nirP^l^ which was in good 
agreement with the terrace structure on bare Ir(lll)P^ 
For intercalation experiments sodium was deposited from 
an outgassed dispenser (SAES getters) keeping the pres- 
sure in the 10~^ Pa regime. 

The plasmons were measured by using a combination 
of high resolution electron loss spectrometer (EELS) as 
electron source and a detector with a low energy elec- 
tron diffraction system (LEED) providing simultaneously 
high energy and momentum (fcy) resolution [28]. Typical 
operating parameters were 25 meV energy resolution at 
a fcy resolution of 1.3 x 10~^A~^. In order to determine 
peak positions and half widths (FWHM), the experimen- 
tal spectra were numerically fitted. Details of the fitting 
procedure can be found in ref.^^Sl. 



III. RESULTS 

As we have already mentioned in recent worlP^ that 
it is not possible to parametrize the observed losses by 
a single loss peak due to their asymmetric form, and it 
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FIG. 1. Experimental set of EEL spectra as a function of 
For monolayer graphene on a) Ir(lll), b) SiC(0001)-6H. 
Solid curves: Fits of the plasmonic losses of the monopole 
(blue) and the multipole plasmons (green). In addition, an 
exponential form of the Drude tail was fitted (not shown). 
For graphene on SiC, also the low energy losses due to Fuchs- 
Kliewer phonons were included in the fits. 



seems that also in such a purely two-dimensional system 
more than one collective electronic excitation is possible. 

This is particularly obvious for graphene on Ir(lll). 
In fig[l] we show a comparison of loss spectra for mono- 
layer graphene grown on Ir(lll) and on the Si face of 
SiC(OOOl), which have been measured for scattering vec- 
tors ranging from zero up to 0.20 A~^. All spectra were 
measured at room temperature at a primary energy of 20 
eV. On the Ir(lll) surface we found apart from asymmet- 
ric peak shapes at low loss energies and small fcy values 
two separable loss peaks that can clearly be discrimi- 
nated at fcy > 0.07 A~^, as seen in fig. ik). The more 
intense peak is assigned to the monopole^heet plasmon. 
Compared to the plasmon of graphene on SiC(OOOl) (see 
fig. [ijj)), the graphene sheet plasmon on Ir(lll) is by a 
factor of five less intense and the full width at half max- 
imum (FWHM) is twice as large (for a detailed analysis 
of the FWHM, see below). 

Nevertheless, the intensities of the second loss peak are 
higher than on SiC(OOOl). While the ratio between the 
sheet plasmon and this high energy mode is 3:1 on SiC 
at fcy = 0.04A-\ it is close to 3:2 on Ir(lll). Quali- 
tatively, and partly even quantitatively, the same results 
were obtained after adsorption of half a monolayer of Na. 
However, the intensity of the sheet plasmon increases 
by roughly a factor of five relative to the elastic peak, 
and the FWHM is reduced. The intensity ratio between 
the two losses again increases to 3:1. This behavior is 
expected for multipole plasmonic excitations when the 
charge distribution normal to the surface varies due to 
different environments^*^. With respect to the systems in- 
vestigated here, the interaction of graphene with Ir(lll) 
is signicantly weaker than for SiC (0001), as directly ob- 
vious from ARPES measurement^^^'^. The Na layer 
turns out to be intercalatecpH. Therefore, it is expected 
not only to change the electronic carrier concentration 
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FIG. 2. a) Plasmon dispersion of the two plasmonic branches of graphene on Ir(lll). b) same type of data evaluation for 
graphene on SiC(OOOl) (n-doped layer n ~ 10^^ cm^^ c) Results for Na intercalated graphene on Ir(lll). For comparison, we 
drew the average linear slopes of a) as dashed lines in the other graphs. 



of graphene, but also to further reduce the interaction 
between graphene and Ir(lll). Therefore, we assign the 
higher loss peak to a multiplasmon mode. In fact, this 
niultipole mode in a 2D sheet may be considered as the 
limiting case of surface plasmon s loca ted at two interfaces 
and interacting with each othei'221^. 

These findings motivated us to further analyze our 
ELS-LEED data of graphene on SiC(OOOl) and to fit the 
data with the same routine and the same procedures as 
for the data on Ir(lll), and to describe the asymmetric 
shape observed there by a second Gaussian peak plus a 
residual background. The results for some of the EELS 
spectra of this system are plotted in fig. [Tja). It is quite 
obvious also in this case that a single loss peak is insuf- 
ficient to describe the shape of the observed loss peaks 
reasonably well. In addition, there is a systematic in- 
crease of separation between the two fitted peaks as a 
function of fcy in both systems. 

Even more informative are plots of the peak positions 
of the losses as a function of /c|| , as shown in fig. [2] In 
fig. [2^), the plot of plasmon dispersion of graphene on 
Ir(lll), we see that the dispersion of both branches is 
strictly linear within error bars. Therefore, also for the 
multipole mode the excitation energy seems to vanish 
in the limit of long wavelengths: w — > for fc|| — > 0. 
This identifies also this second mode as a truely two- 
dimensional mode with a vanishing restoring force in the 
limit of fc|| — ?► 0. 

The ratio between the slopes, i.e. between group ve- 
locities of both branches is 2 ± 0.05. When we now 
turn to the graphene system on SiC(OOOl), again both 
branches are detected (see fig. ^p) , but now the energy 
ratio between these branches is 1.4 ± 0.1. Interestingly, 
the small dip at /c|| — 0.08A~^ exists in both branches. A 
similar evaluation on the Na intercalated Ir(lll) surface 
yields an energy ratio of 1.6 ± 0.07. Qualitatively, and 
even semi-quantitatively, the dispersion of the multipole 
branch follows closely that of the monopole mode. 

Therefore, we will shortly discuss the monopole mode 
here. A more extended description and discussion can 
be found in refs.'^^^. The most remarkable result con- 



cerning the monopole plasmon dispersion is its constancy 
of slope, when substrate and carrier concentration in the 
conduction band are changed. The average slope changes 
by 15% at most, although the doping concentration is 
altered by at least three orders of magnitude from es- 
sentially non-doped in graphene on Ir(lll) to strongly 
doped in the Na intercalated systerrP^. The doping level 
on SiC(OOOl) is in between these two extreme^^. This 
result clear ly co ntradicts the standard theories of a 2D 
electron gaa^'^Il^ which predict a significant dependence 
on carrier concentration, but also on effective masses and 
on the environment. It would be highly accidental, there- 
fore, if these parameters would compensate each other so 
perfectly, in particular for graphene on Ir(lll) with and 
without Na intercalation. It seems that in this gap-less 
electronic system an effective integration over occupied 
electronic states over a range of roughly 1. 5 eV takes 
place yielding an essentially constant effective density of 
electrons participating in plasmonic excitations'^^, pos- 
sibly caused by the short plasmonic life times in these 
systems. However, while experimental parameters like 
step densities on the substrate etc. may modify details 
of the dispersion curves or half widths of the observed 
losses, they only act as second order corrections with this 
effective density still essentially fixed. 

One of these parameters is the exact linearization of 
the dispersion curves observed on Ir, but not on SiC, 
which can be understood by the screening of the plas- 
monic excitation by the bulk metaP^, and has been ob- 
served for 2D plasmons induced by partially filled surface 
states on Be(000lj^and Cu(lll)4^ 

As mentioned, the dispersion of the multipole mode 
shows essentially the same k-dependence as the monopole 
mode. This close relationship may be understood con- 
sidering the restricted possibilities of collective displace- 
ments within a 2D electron gas parallel and normal to the 
plasmon sheet. Especially electronic collective displace- 
ments normal to the graphene sheet, however, should be 
susceptible to the coupling between graphene layer and 
the substrate in addition to the shielding by the sub- 
strate, since coupling to the substrate modifies the restor- 
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FIG. 3. Half widths (FWHM) of loss peaks as a function of k\\, a) for graphene on Ir(lll), b) for graphene on SiC(OOOl), c) 
for Na intercalated graphene on Ir(lll). Squares denote data for the monopole mode, circles those for the multipole mode. 



ing forces. In agreement with this expectation, we find 
that the energetic ratio between monopole and multipole 
mode varies with substrate and in particular with the 
strength of interaction. Considering these ratios for the 
three systems investigated here, however, they obviously 
depend further on the details of interaction and spatial 
electron distribution, not just on the strength of interac- 
tion. On the other hand, the sensitivity of the energetic 
ratio to these parameters is a further indication for the 
existence of the multipole plasmon. 

Coming back to the quantitative properties of the mul- 
tipole plasmon dispersion, we observed even quantitative 
simultaneous changes of both dispersion curves as a func- 
tion of , e.g. the appearance of the dip in both plas- 
mon dispersion curves on SiC(OOOl), can be understood 
as being of the same origin. In this case, the dip was in- 
terpreted as a coupling resonance between plasmonic and 
single particle excitations due to interband transitions^. 
Considering the band structure of graphene, the cross- 
ing of the multipole plasmon dispersion curve and the 
interband transition regime is only marginally shifted to 
smaller fcy values compared to the monopole mode, in 
agreement with the results shown in fig. ^p) and with 
the corresponding resonant increase of the FWHM [see 
fig#)]. 

Generally, we find that the FWHM of monopole and 
multipole plasmons not only exhibit the same functional 
dependence as a function of fcy , they even agree within 
error bar with each other. This is remarkable since with 
their significantly higher excitation energy, the phase 
space for electron hole decay is much larger for the mul- 
tipole plasmons and shorter life times than for monopole 
plasmons are expected. This is obviously not the case and 
indicates that the same decay mechanisms must be effec- 
tive for both modes, which, however, must be of other 
origin. 

An alternative comes from our analysis of the 



monopole plasmon loss and its coupling with struc- 
ture on Ir(lll]p2]_ Here we found that the life time of 
the monopole plasmon on the Ir surface is fully deter- 
mined by the modulation of the graphene layer due to 
its interaction with the substrate. In other words, al- 
though graphene is less strongly coupled to the Ir surface 
than to SiC, so that graphene overgrows the substrate 
stepRoccaOTs in a carpet-like mode, this modulation is 
still sufficient to act as efficient source of momentum 
transfer. Only momentum transfer, however, is neces- 
sary for the transformation of the plasmonic excitation 
to an electron-hole excitation within graphene. Both the 
modulation by the Moire and the step density contribute 
her^2li. From these findings we conclude that this mecha- 
nism with essentially the same probability is also effective 
for the multipole plasmons. 

A similar situation was found for the Na intercalated 
system and long plasmonic wavelengths. Only in this sys- 
tems some deviations between the FWHMs of monopole 
and multipole modes at the shortest wavelengths (/cy > 
were found. For the monopole mode, we ar- 
gued irPS! that the reduced moire-like modulation of the 
graphene film in presence of Na intercalation increases 
the life plasmonic life times compared to the situation 
without Na. This effect is obviously less pronounced for 
the multipole mode, although the life times are still sig- 
nificantly longer than for graphene on the pure Ir(lll) 
surface. 



IV. CONCLUSIONS 

We have demonstrated that in a strictly two- 
dimensional electron system several plasmon modes can 
be excited, and have identified by energy and momentum 
resolved electron loss spectroscopy two of these modes 
for adsorbed single graphene sheets. These were called 
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monopole and multipole plasmonic excitations. The ra- 
tio of excitation energies at constant wavelength seems to 
depend on the type of substrate (metal or semiconduc- 
tor/insulator) and on the coupling strength. More effec- 
tive screening of the excitation by a metallic substrate 
seems to increase the frequency of the multipole plas- 
mon. Reduced coupling between substrate and graphene, 
on the other hand, decreases this ratio. Since reduced 
binding to a substrate is usually coupled with a change 
in average bond lengths, evidenced here by an effectively 
reduced modulation of the graphene layer in the Na inter- 
calated Ir system, the electronic spillout at the interface 
into vacuum is modified. Multipole modes are sensitive 
to the form of the decay of the electron density normal 
to the surface. This explains the observed change of loss 
ratios when the coupling to the substrate is modified. 

Regarding the ratios of intensities, there should be 
clearly a dependence on the oscillator strength of the 
various modes and on their coupling to the multipole 
field of the incident electrons in an EELS experiment. In 
addition, however, there is also a strong modulation of 
intensities by the dynamic electron backscattering proba- 
bility, since this probability varies strongly with incident 
or backscattered electron energy and with their angle of 
incidence, particularly at low energies below 20 eV. As 
an extreme, the observability of certain modes may be 
suppressed at given primary electron energies. 

The overall insensitivity of the measured dispersion 



curves to the actual doping level of graphene is not under- 
stood at present. It clearly invalidates the simple picture 
that only the electron density of the conduction bands 
participates in the plasmonic excitation. Although dop- 
ing induces modifications of the band structure close to 
the Dirac point, which may in part be responsible for 
compensational effects, our experiments show that the ef- 
fective electron density is essentially independent of dop- 
ing concentration. 

The extreme confinement of the electronic excitations 
along a two-dimensional sheet turns out to make these 
excitations extremely sensitive to all kinds of distortions 
of their 2D periodicity. Among these distortions are vari- 
ations of coupling to the substrate by steps and the asso- 
ciated 3D modulation of the graphene film. Our results 
suggest that even for the most perfect substrate of SiC 
with average terrace lengths of more than 100 nm there 
is still a dominance of substrate imperfections for plas- 
mon decay. Therefore, there is an urgent need for even 
more perfect substrates in order to be able to study 2D 
plasmon physics in more detail. 
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